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I.  Introduction 

The  objective  of  this  work  is  to  develop  a  high-order  hybrid  finite  element  technique  that  is 
capable  of  simulating  large  array  antennas  embedded  in  inhomogeneous  media.  Modem  air  and 
land  combat  vehicles  and  warships  are  usually  equipped  with  array  antennas  for  a  variety  of 
communication,  detection,  tracking,  and  surveillance  purposes.  To  reduce  the  radar  signature  and 
any  adverse  effect  on  the  aerodynamic  design,  the  antennas  are  often  conformal  to  the  surface  of 
the  platform  and  sometimes  embedded  in  layered  dielectric  media.  Placing  these  antennas  on  a 
platform  inevitably  introduces  distortion  in  their  radiation  patterns  and  causes  mutual  coupling. 
The  distortion  in  the  radiation  patterns  may  reduce  the  desired  coverage  for  effective 
communications  and  compromise  the  accuracy  for  isolating  and  locating  targets.  The  existence  of 
mutual  coupling,  caused  by  space  waves,  surface  waves,  and  scattering  by  the  platform,  reduces 
the  electromagnetic  isolation  between  the  array  elements  and  consequently  makes  it  difficult  to 
operate  the  array  antennas  simultaneously.  Therefore,  it  is  important  to  develop  accurate 
numerical  prediction  tools  to  characterize  the  radiation  patterns  and  mutual  coupling  of  array 
antennas  mounted  on  a  complex,  often  large,  platform.  As  a  first  step,  it  is  necessary  to  develop 
efficient  tools  that  can  be  used  to  simulate  and  optimize  the  design  of  conformal  array  antennas, 
and  this  is  the  objective  of  the  proposed  research. 


II.  Performed  Research 

Because  of  its  capability  of  modeling  complex  structures  and  inhomogeneous  materials,  the  finite 
element  method  (FEM)  is  best  suited  to  simulating  complex  array  antennas  embedded  in  dielectric 
media  During  the  past  four  years  and  under  the  support  of  this  project,  we  have  developed  a  suite 
of  simulation  tools  based  on  the  FEM  to  deal  with  a  variety  of  array  antennas.  These  include  (1) 
infinitely  large  periodic  phased  arrays,  (ii)  array  antennas  that  are  finite  in  one  dimension  and 
infinitely  periodic  in  the  other  dimension,  (iii)  finite  array  antennas  with  arbitrary  array  elements, 
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and  (iv)  conformal  array  antennas  mounted  on  a  complex,  often  large,  platform.  These  are 
described  in  more  details  below. 

A.  Infinitely  large  periodic  phased  arrays.  The  analysis  of  large  finite  array  antennas  is 
usually  very  time-consuming  and  costly.  Fortunately,  the  characteristics  of  an  array  element  in  a 
large  array  are  similar  to  that  in  an  infinite  array.  Hence,  our  first  step  was  to  develop  a  numerical 
technique  that  could  analyze  infinitely  large  periodic  phased  arrays  embedded  in  a  dielectric 
medium.  In  this  case,  the  Floquet  theorem  is  employed  to  reduce  the  domain  of  analysis  to  a 
single  unit  cell.  Periodic  boundary  and  radiation  conditions  are  enforced  on  the  surface  of  a  single 
unit  cell.  The  asymptotic  waveform  evaluation  (AWE)  technique  is  combined  with  the  FEM  to 
perform  fast  frequency  and  angular  sweeps.  The  computed  antenna  parameters  are  compared  with 
previously  published  results  and  good  agreement  is  obtained.  One  example  concerns  with  a 
circular  patch  array  fed  with  coaxial  lines.  The  periodic  lengths  are  given  by  T_x  =  34  mm  and 
T  y  =  36. 1  mm,  and  the  patch  radius  is  a  =  14.29  mm.  The  outer  and  inner  radii  of  the  coaxial  line 
used  in  the  simulation  are  1.492  and  0.456  mm,  respectively.  The  Sn  parameter  at  the  coaxial  port 
is  shown  in  Figure  1.  Both  the  magnitude  and  phase  of  the  calculated  Sn  parameter  agree  with  the 
results  obtained  by  the  hybrid  generalized  scattering  matrix  (GSM)  and  the  FEM. 


Figure  T.  Magnitude  and  phase  of  the  Sn  parameter  for  the  microstrip  patch  array  consisting  of 
circular  patches  with  radius  a  =  14.29  mm,  substrate  thickness  h  =  0.79  mm,  and  relative 
permittivity  of  2.33.  The  feed  is  modeled  using  a  precise  coaxial  line  model. 

B.  Finite-by-infinite  array  antennas.  In  this  case,  the  array  is  infinite  and  periodic  in  one 
dimension  and  finite  and  arbitrary  in  the  other  dimension.  The  analysis  of  this  type  of  array  is 
important  because  it  enables  fast  and  efficient  investigation  of  the  truncation  (or  edge)  effect  of  a 
finite  array,  which  distinguishes  a  finite  array  from  an  infinite  array.  Here,  the  Floquet  theorem  is 
employed  to  reduce  the  domain  of  analysis  to  one  row  of  array  elements.  Again,  periodic 
boundary  and  radiation  conditions  are  formulated  for  a  unique  FEM  solution  of  the  fields  inside 
the  domain  of  analysis.  As  an  example,  consider  a  5  x  infinite  vivaldi  array  antenna.  Figure  2 
shows  the  unit  cell  geometry  of  the  antenna  and  the  normalized  radiation  power  pattern  in  the  E- 
plane  as  a  function  of  the  receiving  angle  at  3, 4,  and  5  GHz. 
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Figure  2:  Geometry  of  the  vivaldi  array  antenna  (a  single  element)  and  the  E-plane  radiation 
pattern  of  the  5  x  infinite  vivaldi  array  antenna. 


C.  Arbitrary  finite  array  antennas.  Among  the  three  types  of  arrays  considered  here,  this  is 
most  challenging  because  even  for  an  array  antenna  of  moderate  size,  the  computational 
requirements  become  very  excessive.  This  is  particularly  true  when  the  computational  domain  is 
truncated  using  a  boundary  integral  equation,  which  yields  a  fully  populated  submatrix  for  the 
fields  over  the  entire  aperture  of  a  finite  array.  To  alleviate  this  problem,  the  adaptive  integral 
method  (AIM)  is  employed  to  efficiently  evaluate  the  boundary  integrals  with  the  aid  of  the  fast 
Fourier  transform  (FFT).  Figure  3  shows  the  calculation  of  the  mutual  coupling  (S-parameters) 
among  a  2  x  2  microstrip  patch  array  antenna.  More  specifically,  it  gives  the  plots  of  the  4  x  4 
scattering  matrix  from  1  to  3  GHz.  The  FEM  results  obtained  with  the  boundary  integral 
truncation  (circles)  compared  very  well  with  the  results  (lines)  of  another  calculation  by  the  time- 
domain  FEM  using  perfectly  matched  layers  for  mesh  truncation. 


D.  Conformal  antennas  mounted  on  a  large,  complex  platform.  All  the  antennas  are 
eventually  mounted  on  a  platform  such  as  air  and  sea-going  vehicles.  The  interaction  between  the 
antennas  and  platform  is  a  critical  factor  in  the  antenna  design;  unfortunately,  it  also  represents 
one  of  the  most  challenging  problems  faced  by  the  computational  electromagnetics  community.  In 
this  project,  we  have  developed  a  hybrid  technique  that  couples  the  finite  element  method  and  a 
surface  integral  equation  for  evaluating  conformal  antennas  on  a  large,  complex  platform.  High- 
order  curvilinear  elements  were  employed  to  accurately  represent  the  electric  and  magnetic  fields 
as  well  as  the  antenna  geometry.  The  high  computational  complexity  associated  with  the  solution 
of  the  surface  integral  equation  was  alleviated  by  the  use  of  the  multilevel  fast  multipole  algorithm 
(MLFMA).  A  physical-based  preconditioner  was  designed  to  drastically  decrease  the  number  of 
iterations  required  for  convergence.  Figure  4  shows  a  conformal  microstrip  patch  antenna  housed 
in  a  cavity  that  resides  on  a  platform  consisting  of  a  conducting  circular  cylinder  and  a  conducting 
plate  (wing).  Also  shown  are  the  measured  data  by  Mission  Research  Corporation,  and  it  is  seen 
that  the  numerical  results  agree  well  with  the  measurement  for  both  co-  and  cross-polarizations. 
Using  the  physical-based  preconditioner,  the  final  solution  was  obtained  with  only  10  iterations. 
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This  example  clearly  demonstrates  the  significant  effect  of  the  platform  on  the  antenna’s  radiation 
patterns. 


Figure  3:  Scattering  matrix  S(f)  of  the  2x2  microstrip  patch  array  antenna.  The  graph  on  the  i-th 
row  and  j-th  column  shows  ISjj(f)l  for  the  frequency  band  from  f  =  1  to  3  GHz. 


Figure  4:  (a)  Cavity-backed  patch  antenna  mounted  on  a  cylinder  with  a  wing,  (b)  H-plane 
radiation  pattern  at  3.3  GHz.  The  simulation  is  performed  using  the  higher-order,  hybrid 
FEM/MLFMA  method. 
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The  numerical  techniques  developed  in  this  project  all  have  the  following  important 
characteristics:  (i)  higher-order  geometrical  modeling,  (ii)  higher-order  field  discretization,  (m) 
hybridization  with  surface  integral  equations  using  fast  algorithms,  (iv)  a  highly  effective 
preconditioner,  and  (v)  accurate  antenna  feed  modeling.  These  are  described  in  more  details 

below. 


E.  Higher-order  geometrical  modeling.  Higher-order  geometrical  modeling  is  important  in 
order  to  model  curved  surfaces  accurately  so  that  surface  waves  can  be  simulated  with  a  high 
accuracy.  Fortunately,  higher-order  geometrical  modeling  can  be  achieved  rather  easily  using 
parametric  mapping,  which  can  transform  a  planar  triangular  patch  into  a  quadratic  or  cubic 
curved  triangle.  The  same  holds  for  tetrahedral  elements.  This  allows  the  field  expansion  and 
discretization  to  be  carried  out  on  a  regular  planar  triangular  patch  or  inside  a  regular  tetrahedron 

in  the  transformed  space. 

F.  Higher-order  field  discretization.  Higher-order  field  discretization  is  critical  for  achieving 
a  higher-order  convergence;  or  in  other  words,  to  obtain  a  highly  accurate  solution  wit  a 
minimum  number  of  unknowns,  which  translates  to  a  better  efficiency.  Higher-order  fie  d 
discretization  is  particularly  important  for  simulating  wave  phenomena  because  it  can  greatly 
suppress  the  accumulative  grid  dispersion  error.  In  this  work,  we  employed  bo5*  ^er-or  er 
interpolatory  and  hierarchical  vector  basis  functions  in  the  implementation  of  the  FEM.  We  have 
performed  a  systematic  convergence  study,  which  confirmed  higher-order  convergence  as 
expected.  The  use  of  hierarchical  basis  functions  permits  the  future  application  of  hp-adaptive 
refinement  schemes  that  would  yield  optimal  convergence  for  specified  accuracy. 


G  Hybridization  with  surface  integral  equations  using  fast  algorithm.  As  is  well  knovm, 
the  most  challenging  problem  for  the  FEM  to  model  unbounded  electromagnetic  problems  u;  the 
mesh  truncation.  Although  the  use  of  absorbing  boundary  conditions  (ABC),  including  perfectly 
matched  layers  (PML),  is  straightforward,  the  resulting  solution  becomes  approximate  because  no 
ABC  and  PML  can  simulate  a  perfect  absorbing  boundary  (which  has  no  non-physical  reflection). 
The  best  approach  is  to  use  surface  integral  equations  that  are  exact  and  can  be  placed  ect  y  on 
the  surface  of  the  object  to  minimize  the  computational  domain  (thus  the  accumulative  grid 
dispersion  error).  However,  the  use  of  a  surface  integral  equation  on  the  truncation  surface  results 
in  a  full  submatrix,  whose  computation  becomes  very  expensive.  If  we  let  N  to  denote  the  number 
of  unknowns  associated  with  the  surface  integral  equation,  the  computational  complexity  of  the 
submatrix  is  proportional  to  0(N2).  This  has  long  been  regarded  as  the  bottleneck  of  the  surface 
integral  equation.  In  this  work,  we  removed  this  bottleneck  and  reduced  the  computational 
complexity  to  0(N  log  N)  by  using  either  AIM  or  MLFMA,  depending  on  the  specific  type  of 

antennas  to  be  analyzed. 

H.  A  highly  effective  preconditioner.  As  mentioned  earlier,  the  use  of  a  surface  integral 
equation  on  the  truncation  surface  results  in  a  full  submatrix,  which,  when  combined  with  t  e 
FEM  matrix,  yields  a  partly  sparse  and  partly  full  system  matrix.  This  system  matrix  ^rather  ill 
conditioned  and  a  direct  application  of  an  iterative  algorithm,  such  as  those  based  on  the  Kry  ov 
subspace,  converges  very  slowly.  We  have  tested  a  few  algebraic  preconditioners;  but  the 
improvement  was  not  significant  enough  to  yield  an  efficient  solution.  For  this  wmk,  we 
developed  a  novel,  physical-based  preconditioner,  which  turned  out  to  be  very  highly  effective. 
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This  preconditioner  is  based  on  our  discovery  that  the  spectrum  of  a  FEW  ABC  system  is  very 
similar  to  that  of  the  FEM/CFIE,  where  CFIE  stands  for  the  combmed  surface  integral  equation 
which  is  the  most  accurate  surface  integral  equation  used  for  the  mesh  truncation.  As  a  result,  t 
FEM/ABC  system  serves  as  a  very  effective  preconditioner  for  the  FEM/CFIE  matrix  equation. 
With  this,  convergence  is  usually  achieved  within  a  few  iterations.  This  precon  ltioner  app  es 
both  antenna  and  scattering  problems,  and  its  development  is  considered  a  major  progress  in  the 
FEM  analysis  of  unbounded  electromagnetic  problems. 

/  Accurate  antenna  feed  modeling.  One  of  the  very  critical  aspects  in  the  development  of 
any  numerical  methods  for  accurate  antenna  analysis  is  the  modeling  of  ante^^,^°  | 
reasonably  accurate  results  for  radiation  patterns  can  be  obtained  using  simjMed  feed  m 
accurate  characterization  of  input  impedance  and  mutual  coupling  (mutual  impedances  or  S 
parameters)  can  only  be  obtained  using  accurate  feed  models.  In  practice,  most  antennas  are  fed 
by  coaxial  lines  However,  most  numerical  analyses  use  simplified  models  to  make  the  ProWem 
more  tractable.  For  example,  in  the  integral  equation  based 

differential  equation  based  FEM,  the  delta  voltage  gap  and  current  probe/filament  are  often  used 
““Tan 2m* .  Although  these  staple  models  axe  capable  of  generating  accurate  radratron 
patterns  and  in  some  cases  input  impedances  for  simple  wire  and  patch  antennas,  they  often  do 
SHSnS  for  complex^  antennas  and,  worst  of  all,  their  modelmg  errors  cannot  be 
systematically  reduced  (controllable).  Although  the  use  of  a  magnetic  frill  model  improves  t 
accuracy  of  Lpedance  calculation  for  wire  antennas,  its  practical  application  is  still  very  limited 
In  ^  work,  we  have  developed  an  accurate  method  to  model  a  variety  of  antenna  feed.  This 
model  is  directly  based  on  the  full-wave  analysis  of  electric  and  magnetic  fields  in  the  fee 
^ures,  Stead  of  using  the  voltage  and  current  concepts.  This  model  involves  a  reference 
surface  often  chosen  to  be  close  to  the  coaxial  or  waveguide  opening,  and  then  represents  t  e 
total  field  as  the  superposition  of  the  incident  and  reflected  waves  (mcludmg  h*i-o* ^odes 
excited  by  the  structures  in  the  vicinity  of  the  feed).  An  exact  boundary  condition  can  then  be 
derived  which  can  be  incorporated  into  the  FEM  solution  of  the  antenna  problem.  As  a  result  the 
^  for  the  numerical  antdysis  is  the  incident  field  in  the  fad  line  and  the  output  rs  the  reflection 

coefficient  (S„),  from  which  the  input  impedance  can  be  readdy  ,C“°  e J”  *t“^n be 
■....triple  antennas,  which  is  of  particular  interest  in  this  project,  a  full  S-parameter  matra  can  te 
computed  by  exciting  each  antenna.  This  model  is  applicable  to  most  antenna  feeds  (coaxial  hn  , 
waveguidesfand  other  transmission  lines).  Except  for  numerical  discretization,  the model Its exact 
and  its  numerical  solution  is  error  controllable.  More  important  since  dns .wM plates 
reflected  waves  in  the  feed  lines,  one  can  incorporate  the  effect  of  feed  structures,  such 
corporate  feed,  into  the  array  antenna  analysis.  All  the  results  shown  above  were  obtained  using 

this  feed  model. 


m.  Beyond  the  Proposed 

In  addition  to  the  work  described  above,  which  has  successfully  fulfilled  the  proposed  objective 
we  have  started  to  explore  two  topics  that  are  closely  related  to  the  proposed  project  CO  As 
mentioned  earlier,  an  efficient  and  accurate  mesh  truncation  technique  is  critical  to  the  FEM 
analysis  of  large-scale  antenna  problems.  In  this  work,  we  have  employed  surface  integral 
equations  with  fast  algorithms  for  this  purpose.  Researchers  in  computational  acoustics,  however, 
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reported  better  performance  of  infinite  elements  over  surface  integral  equations.  This  has  also  led 
some  researchers  in  computational  electromagnetics  to  pursue  the  development  of  infinite 
elements  for  the  mesh  truncation.  We  have  developed  infinite  elements  and  conducted  a  detailed 
comparative  study  against  the  FEM  with  a  surface  integral  equation  and  concluded  that  the 
second  approach  is  more  accurate  and  efficient  than  the  first  one  using  infinite  elements.  We  also 
discovered  that  for  three-dimensional  electromagnetic  problems,  infinite  elements  can  only  be 
developed  for  a  spherical  surface  because  of  the  vectorial  nature  of  the  field  wave  equation, 
although  for  the  acoustic  field  infinite  elements  have  been  developed  for  spheroidal  surfaces.  The 
use  of  spherical  truncation  surfaces  is  highly  inefficient  and  has  little  practical  use  for  large-scale 
antenna  problems,  whereas  surface  integral  equations  can  be  placed  right  on  the  surface  of  the 
antenna  to  be  simulated.  (2)  Recognizing  the  importance  and  urgent  need  to  develop 
computational  techniques  to  simulate  ultra-wide  band  antennas,  we  started  to  develop  the  higher- 
order  time-domain  finite  element  method  (TDFEM).  Beside  its  capability  to  perform  a  broadband 
analysis,  the  TDFEM  can  effectively  model  nonlinear  and  inhomogeneous  devices/components/ 
medium.  Its  geometrical  modeling  capability  is  excellent  and  the  time-stepping  can  be  made 
unconditionally  stable.  Hence,  this  method  overcomes  all  the  shortcomings  of  the  frequency- 
domain  and  FDTD  techniques.  Its  major  drawbacks  are  (i)  the  need  to  solve  a  large  sparse  matrix 
equation  in  each  time  step  and  (ii)  the  implementation  of  perfectly  matched  layers  (PML)  becomes 
very  complicated.  However,  the  development  of  large,  sparse  matrix  solvers  has  partially 
alleviated  the  first  drawback.  We  have  concentrated  our  effort  to  develop  a  stable  PML 
formulation  for  the  mesh  truncation  of  the  TDFEM.  We  are  convinced  that  the  TDFEM  offers  a 
promising  approach  to  the  simulation  of  ultra  wide  band  antennas  with  nonlinear  active  devices. 
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Technical  Significance: 

□  Fast  and  accurate  numerical 
simulation  is  critical  to  the  design  of 
advanced  complex  array  antennas 
embedded  in  dielectric  media 

□  Characterization  of  the  mutual 
coupling  is  required  for  the  design  of 
multifunction  (transmitter/receiver) 
antenna  systems 

Impact: 

□  Simulation  of  large,  complex  antenna  systems  for  navy  and  other  applications 

□  Advancement  of  knowledge  in  computational  electromagnetics  for  complex 
geometries 

□  Development  of  new  techniques  and  algorithms  that  hybridize  different 
numerical  methods 

□  Methodologies  also  applicable  to  scattering,  wave  interaction,  and  circuits 
problems 
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Introduction 


A  complete  set  of  FEM-based  numerical  simulation 
techniques  has  been  developed  for  a  systematic  analysis  of. 

■  Infinite  periodic  phased  arrays 

■  Single  unit  cell  analysis  for  basic  characteristics 

■  Finite-by-infinite  array  antennas 

■  Study  of  truncation  (edge)  effects 

■  Arbitrary  finite  array  antennas 

■  Highly  challenging  for  a  large  finite  array 

■  Conformal  antennas  mounted  on  a  platform 

■  The  ultimate  goal  -  Antenna/platform  interaction  modeling 
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Common  characteristics  of  the  FEM-based  numerical  simulation 
techniques  developed  here: 

■  Higher-order  geometrical  modeling 

■  Higher-order  field  discretization 

■  Accurate  antenna  feed  modeling 

■  Boundary  integral  equations  (BIE)  as  exact 
mesh  truncation 

■  Fast  algorithms  (FMM  &  AIM)  for  BIE  evaluation 

■  Effective  preconditioner  to  accelerate 
convergence  of  iterative  solution 
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Feed  Modeling 


Accurate  feed  modeling  is  critical 
for  antenna  simulation  (especially 
for  characterization  of  input 
impedance,  mutual  coupling,  &  S 
parameter). 

Feeding  structure  included  by 
coaxial  cables  with  waveguide 
ports  (rigorous  formulation). 


Full-wave  computation. 
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Infinite  Periodic  Array 
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□  Why  investigate  finite  by  infinite  (FBI)  arrays? 

-  A  good  model  to  study  the  effects  of  array  truncation 

—  Numerically  less  demanding  compared  to  finite  by  finite  arrays 

-  Approximation  to  some  realistic  applications 


Dielectric 


PEC 


Spectral  domain  formulation: 

H— (r)  =  zA^L  JV**'  JG(*V  ,*,)•  M(kxp ,  ky )  eA’'dk) 


PEC\ 


Conductor 


Mixed  spectral-spatial  domain 
formulation: 

Hsc°'(r)  =  -2jk0Y0  Yjelt,x \G{kVJ,y\y')mk^y')dy' 
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Fin ite-by-ln  finite  Array 


•  Frequency:  5  GHz 

•  Unit  cell  dimension:  32  X  34  X  33.3  mm 

•  Reflection  calculated  at  individual  coaxial 
ports  for  E-plane  scan 


■  infintle  array 

•  5-by-lnflnlte  array,  column  1 
-  5-by-ln!lnlte  array,  column  2 

■  5-by-lnfinlte  array,  column  3 
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•  Reflection  calculated  at  individual  coaxial 
ports  for  broadside  radiation 


Radiation  Power  Pattern 
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Geometry  of  the  vivaldi  array  antenna  (a 

single  element)  and  the  E-plane  radiation  <*  = 33.3mm 

pattern  of  the  5  x  infinite  vivaldi  array  w=30mm 
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Adaptive  integral  method: 

-  Two  sets  of  meshes 

-  Auxiliary  sources  equivalently 
represent  the  unknowns 

-  Use  of  FFTs 
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Finite  Array  Antennas 


Radiation  pattern  from 
antenna  element  #1 
Normalized  Poynting 
vector  extends  from 
origin  to  the  colored 
surface 
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Finite  Array  Antennas 
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Antenna/Platform  Interaction 


Final  combined  matrix 
system: 

1 .  The  system  is  usually  very  large; 
hence,  iterative  solution  is  the 
only  choice. 

2.  Unfortunately,  convergence  is 
very  slow. 
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ABC-Based  Preconditioner: 


Replace  CFIE  by  1st-order  ABC: 

n0  x  E(r)  +  Z0no  xh0x  H(r)  =  0  r  £  S0 

Combined  system  becomes: 
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Note:  1.  A  large,  purely  sparse  matrix.  sparse  matrices 
2.  A  very  effective  preconditioner. 
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ABC-Based  Preconditioner 
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Note:  When  higher-order  elementsO 
are  used,  the  convergence  without 
preconditioning  is  even  slower, 
whereas  the  convergence  with 
preconditioing  remains  the  same. 
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Conclusion 


A  complete  set  of  FEM-based  simulation  techniques: 

■  Infinite  periodic  phased  arrays 

■  Finite-by-infinite  array  antennas 

■  Arbitrary  finite  array  antennas 

■  Conformal  antennas  mounted  on  a  platform 
Common  characteristics: 

■  Higher-order  geometrical  modeling 

■  Higher-order  field  discretization 

■  Accurate  antenna  feed  modeling 

■  BIE  as  exact  mesh  truncation 

■  FMM  &  AIM  for  BIE  evaluation 

■  Effective  preconditioner 
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